First principles study of the electronic and magnetic structures of the tetragonal and 

orthorhombic phases of Ca3Mn207 
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On the basis of density functional theory electronic band structure calculations using the aug- 
mented spherical wave method, the electronic and magnetic properties of the orthorhombic and 
tetragonal phases of Ca3Mn207 were investigated and the spin exchange interactions of the or- 
thorhombic phase were analyzed. Our calculations show that the magnetic insulating states are more 
stable than the non-magnetic metallic state for both polymorphs of Ca3Mn207, the orthorhombic 
phase is more stable than the tetragonal phase, and the ground state of the orthorhombic phase 
is antiferromagnetic. The total energies calculated for the three spin states of the orthorhombic 
phase of Ca3Mn207 led to estimates of the spin exchange interactions J,j„ = -3.36 meV and J„„„ = 
-0.06 meV. The accuracy of these estimates were tested by calculating the Curie- Weiss temperature 
within the mean-field approximation. 

PACS numbers: 62.25.-|-g 71.15.Mb 71.20.Nr 
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I. INTRODUCTION 



Due to their potential technological application, mag- 
netic manganese oxides have received much attention in 
the past few years [1-3]. It is an important theoret- 
ical issue to understand their electronic and magnetic 
properties [4, 5]. In accounting for such properties of 
magnetic oxides, density functional theory (DFT) [6, 7] 
electronic structure calculations using the local spin den- 
sity approximation (LSDA) [8] and the generalized gradi- 
ent approximation (GGA) [9] have proven very reliable. 
In the Ruddlesden-Popper (RP) family of manganites 
Ca„+iMn„03„+i [10], the diamagnetic CaO rocksalt lay- 
ers alternate with the magnetic (CaMn03)„ layers along 
the crystallographic c-direction, where n represents the 
number of MnOa perovskite sheets in each (CaMnOs),! 
layer. To emphasize this structural feature, the formulas 
for these oxides can be rewritten as (CaO)(CaMn03),i 
with n = 1, 2, 3, ...,00 . Known examples of this RP 
family are Ca2Mn04 (n = 1), Ca3Mn207 (n 2) and 
Ca4Mn30io (n = 3) [11]. In a recent electronic structure 
study [12], it was shown that the n = 1 RP member, 
Ca2Mn04, and its reduced phase Ca2Mn03.5 are anti- 
ferromagnetic (AF) insulators with band gap of approxi- 
mately ~ 1 eV. According to a recent neutron diffraction 
study [13] at room temperature, the n = 2 RP mem- 
ber, Ca3Mn207, exists as a mixture of the majority or- 
thorhombic and the minority tetragonal phases. The 
main structural difference between the two is depicted 
in Figure 1. Due to the tilting of the MnOe octahe- 
dra, the Mn — O — Mn angles between two adjacent - 
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corner sharing- octahedra are considerably smaller 180° 
(i.e., 159.1, 163.2 and 163.9°). Ca3Mn207 undergoes a 
three-dimensional (3D) magnetic ordering below Tn ^ 
110 - 115 K with a G-type AF structure [13]. As a part 
of our continuing studies on calcium manganese oxides 
[4, 12, 14, 15], we examine in the present work the elec- 
tronic and magnetic structures of the two polymorphs of 
Ca3Mn207 on the basis of first principles DFT electronic 
structure calculations. 



II. CALCULATIONS 

Our DFT calculations employed the all-electron aug- 
mented spherical wave (ASW) method in its scalar- 
relativistic implementation[16, 17]. As preliminary cal- 
culations revealed that the GGA describes the effects of 
exchange and correlation better than the LSDA, hence 
leading to a more accurate description of the energy dif- 
ferences between different ordered magnetic states, we 
opted for a use of the GGA in the form proposed by 
Perdew, Burke, and Ernzerhof [9]. In the ASW method, 
the wave function is expanded in atom-centered aug- 
mented spherical waves, which are Hankel functions and 
numerical solutions of Schrodinger's equation, respec- 
tively, outside and inside the so-called augmentation 
spheres. In order to optimize the basis set, additional 
augmented spherical waves were placed at carefully se- 
lected interstitial sites. The choice of these sites as well 
as the augmentation radii were automatically determined 
using the sphere-geometry optimization (SGO) algorithm 
[18]. The Ca 4s, Mn 4s, Mn 3d and O 2p states were 
treated as valence states, and the low-lying O 2s states 
as core states. The basis set was complemented by wave 
functions of s-, p-, and possibly, d-symmetry centered 
at the interstitial sites. A sufficiently large number of k 
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points was used to sample the irreducible wedge of the 
Brillouin zone, and the energy and charge differences be- 
tween successive iterations were converged below AE = 
10~^ Ryd. and A Q = 10~^, respectively, to obtain ac- 
curate values of the magnetic moments and accurate to- 
tal energy differences between various ordered magnetic 
states. To extract more information about the nature 
of the interactions between the atomic constituents from 
electronic structure calculations, the crystal orbital over- 
lap population (COOP) [19] or the crystal orbital Hamil- 
tonian population (COHP) [20] may be employed. Both 
approaches provide a qualitative description of the bond- 
ing, nonbonding and antibonding interactions between 
two atoms. A slight refinement of the COHP was re- 
cently proposed in form of the "energy of covalent bond" 
(ECOV), which combines COHP and COOP to calculate 
quantities independent of the choice of the zero of poten- 
tial [21]. Both COOP and ECOV give similar general 
trends, but COOP, when defined within plane-wave ba- 
sis sets, exaggerates the magnitude of antibonding states. 
In the present work the ECOV was used for the chemical 
bonding analysis. In the plots, negative, positive and zero 
ECOV magnitudes arc relevant to bonding, antibonding 
and nonbonding interactions respectively. 



III. RESULTS AND DISCUSSION 
A. Non-magnetic configuration 

To gain insight into the chemical bonding in the two 

forms of Ca3Mn207, their non-magnetic (NM) states 
were calculated by performing non-spin-polarized elec- 
tronic structure calculations by enforcing spin degener- 
acy for all species. Our calculations show charge trans- 
fer from the cation sites (Ca^"*", Mn^+) to the anionic 
sites of oxygen as well as toward the empty space, but 
the amount of charge transfer is not significant of an 
ionic character of the atomic constituents. The plots of 
the projected density of states (DOS) for the constituent 
atom sites are presented in Figure 2. The projected DOS 
plots calculated for the five Mn 3d orbitals are presented 
in Figure 3, and the ECOV plots calculated for the vari- 
ous Mn-0 bonds in Figure 4. 

The comparison of Figures 2 and 4 show that the va- 
lence bands (VB's) are dominated by Mn-0 bonding. 
The Mn-0 interactions are bonding in the major part 
of the VB's with strong contributions in the -8, -6 eV 
region (Figure 2), i.e., where the Mn e^-orbitals (Figure 
3) make a-bonding with the O 2p orbitals. The s-like Ca 
states are smeared into the whole range of the VB's (Fig- 
ure 2). For both polymorphs of Ca3Mn207, the Fermi 
level Ep occurs near the DOS peak dominated by the 
Mn 3d states so that the DOS value at the Fermi level, 
n(Eir), is large. The DOS peak around Ep arises from 
the t2p-block bands in which the Mn t2p-orbitals make 
weak TT- antibonding interactions with the O 2p orbitals. 
The existence of Mn-0 antibonding in this DOS peak is 



clearly seen from the ECOV plot (Figure 4). The occur- 
rence of the Fermi level around a sharp DOS peak with a 
large n(Ei?) value, indicates an instability towards a fer- 
romagnetic state. This can be inferred from the Stoner's 
mean field theory analysis (see [22] for a review on mag- 
netic oxides) whereby a nonzero magnetic moment on 
the Mn sites should occur when two spin populations are 
accounted for in the calculations. 



B. Magnetic states 

Spin-polarized electronic structure calculations were 
carried out for the ferromagnetic (FM) state of the or- 
thorhombic and tetragonal phases of Ca3Mn207, and for 
the AF state of the orthorhombic phase of Ca3Mn207. 
There are a number of possible AF spin arrangements. 
In the present study, we considered two types of AF ar- 
rangements, i.e., the G- and A-type AF states shown in 
Figure 5. In the G-type AF state the antiferromagneti- 
cally ordered planes of Mn sites arc antiferromagnetically 
coupled. In the A-type AF arrangement, the ferromag- 
netically ordered planes of Mn sites are antiferromagnet- 
ically coupled within each perovskite layer. 

Results of our calculations for the FM state of the 
tetragonal and orthorhombic phases of Ca3Mn2 07 are 
summarized in Table I as well as Figure 6. For both 
phases the FM state is substantially more stable than 
the NM states (i.e., EpM < E^vm), and the total mo- 
ment per formula unit (FU) is 6.0 (Table I), as ex- 
pected for the high-spin Mn^+ (d^) ions inCa3Mn207. 
The total energy obtained for the orthorhombic phase is 
by 170 meV/FU lower than that of the tetragonal phase, 
which is consistent with the finding that the orthorhom- 
bic and tetragonal phases are the majority and minority 
phases, respectively. The projected DOS plots calculated 
for the FM state of the orthorhombic Ca3Mn207 are pre- 
sented in Figure 6a. The up-spin and down-spin bands 
are strongly split such that the up-spin t2c,-block bands 
are completely filled while the down-spin t2g-block bands 
are empty, as expected for the high-spin Mn^+ (d'') ions. 
The down- and up- spin bands have a band gap of about 
1.0 and 0.3 cV, respectively. Given the fact that both 
the LSDA and the GGA underestimate band gaps, the 
FM state of Ca3Mn207 should be insulating experimen- 
tally. However, the FM state is not the magnetic ground 
state of Ca3Mn2 07, as will be discussed in next section. 
Figure 6b shows the ECOV plots for the Mn-0 bonds 
calculated for the up-spin and down-spin bands of the 
orthorhombic phase. The VB's exhibit Mn-0 bonding 
interactions up to -2 eV in both up-spin and down-spin 
bands. The energy region beteween -2 eV and the Fermi 
level shows Mn-0 antibonding interactions due to the 
TT-type antibonding interactions of the t2g-block bands 
in the up-spin bands, but this feature is missing in the 
down-spin bands because the down-spin t2g-block bands 
are raised above the Fermi level. Thus, the exchange 
splitting in the FM state enhances Mn-0 bonding inter- 
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actions in Ca3Mn2 07. 

Results of calculations for the G- and A- type AF state- 
ments of the orthorhombic phase of Ca3Mn207 are sum- 
marized in Tabic II and Figure 7. Both AF states arc 
more stable than the FM state, and the G-type AF state 
is slightly more stable than the A-type AF state (Ta- 
ble II). The latter is consistent with the neutron diffrac- 
tion study [13], which found a G-type AF ordering for 
Ca3Mn207. While the resulting magnetization is zero as 
expected in AF order, the atomic magnetic moments and 
the resulting t,J, total magnetization are reduced with 
respect to their value in the ferromagnetic configuration. 
This is due to the symmetry breaking when the AF lat- 
tice was accounted for. The DOS and band structure are 
shown in Figure 7 for one magnetic sublattice. The simi- 
larities with the features observed in Fig. 6 for the ferro- 
magnetic case are present here too. So at least from the 
point of view of the electronic band structure some rela- 
tionship is persistent although the AF state is the ground 
state with an insulating character. From the band struc- 
ture plot (Figure 7b) the large dispersion of the bands 
within the CB illustrates furthermore the itinerant char- 
acter of the Cg states. The gap of ^ 0.4 eV is found 
direct, between Ty and Tc- Since both the LSDA and 
the GGA underestimate the band gap, the system can 
be expected to be insulating experimentally. 



E{FM) - E{G - AF) = ■il.Sm.eV/Mn, 
E{A - AF) - E{G - AF) = 29.7meV/Mn. (2) 

From Eqs. (1) and (2), it is estimated that J„„ = -3.36 
meV and J„„„ = -0.06 meV, i.e., = -39 K and ^^f^ = 
-0.7 K, where is the Boltzmann constant. The in- 
teraction is strongly AF, and the NNN interaction is very 
weakly AF. Therefore, within the mean-field approxima- 
tion, the Curie- Weiss temperature of Ca3Mn207 can be 
estimated from the J„„ value according to the expression 
(3) 

_ zJnnSjS + 1) , , 

where z = 5, and S = | for the high-spin Mn^+. For 

Jnri/k_B ^ -39 K, we obtain 8= -244 K, which is some- 
what smaller in magnitude than the experimental value 
©exp = - 465 K [26]. This suggests that our calcula- 
tions underestimated the spin exchange interactions of 
Ca3Mn207 by a factor of approximately two. However 
such a magnitude is expected in the framework of the 
approximation used within the DFT and its functionals 
when comparisons with mean-field analysis results are 
carried out. 



IV. ANALYSIS OF THE SPIN EXCHANGE 
INTERACTIONS 



It is of interest to estimate the nearest-neighbor (NN) 
spin exchange J„„ and the next-nearest-neighbor (NNN) 
spin exchange J„„„ interactions of the (CaMn03)2 
double-perovskite layer in the orthorhombic phase of 
Ca3Mn207. The spin exchange interactions between the 
adjacent double-perovskite layers are expected to be neg- 
ligible, so we consider only those interactions within a 
double-perovskite layer. Then, each Mn site has five NN 
and eight NNN spin exchange interactions. For simplic- 
ity, we assume that all five NN interactions are identical, 
and so are all the eight NNN interactions. Then, in terms 
of inn and J„„„, the energies per Mn site of the FM, G- 
type AF and A-type- AF states are written as [23-25] 



E{FM) = - (-5J„ 



iV2 ^2 

8J„„„)( — ) = ( — )(-2.5J„ 



V. CONCLUDING REMARKS AND OUTLOOK 

Our electronic structure study reveals that for both 

polymorphs of Ca3Mn207, the magnetic states are more 
stable than the non-magnetic states. The spin polar- 
ization makes Ca3Mn207 insulating and lowers its en- 
ergy mainly through ex(;hange leading to a magnetic mo- 
ment on Mn expected for a high-spin Mn^+ (d"^) ion. 
The down-spin states provide more Mn-0 bonding than 
do the up-spin states because the down-spin Mn t2g- 
block bands are raised above the Fermi level. In agree- 
ment with experiment, our calculations show that the 
orthorhombic phase is more stable than the tetragonal 
phase, and the ground state of the orthorhombic phase 
is G-type AF. The total energies of the FM, A-type AF 
and G-type AF states of Ca3Mn207 obtained from our 
calculations lead to the estimates of the spin exchange 



mteractions 



-39 K and 



-0.7 K. The Curie- 



E{G -AF) = - (5 J„„ - 8 J, 



nnn)(. . ) — ( . )^-^Jnn 
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Weiss temperature Q estimated fi'om these spin exchange 
parameters is approximately half the experimental value, 
4 JnSeX^^at our calculations underestimated the values of the 
spin exchange parameters by a factor of approximately 

A T twp. 



E{A-AF) = i (-3J„„)(^) = (^)(-1.5J„„),(l) 
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TABLE I: Total energies and magnetic moments of the atomic 
sites calculated for the ferromagnetic states of the tetragonal 
and orthorhombic phases of Ca3Mn2 07. The energies are 
given in units of eV per half the FU (i.e., per Mn) with re- 
spect to those of the non-spin-polarized state. The magnetic 
moments are given in units of Bohr magneton (/is). 



Parameters 



Tetragonal Orthorhombic 



AEfai-nm 


(eV) 


- 1.442 


-1.610 


AEfai-nm 


(eV)/Mn 


- 0.721 


-0.831 


^(Mn) 




2.729 


2.719 


m(oi) 




0.088 


0.092 


m(02) 




-0.001 


0.047 


m(03) 




0.073 


0.056 


m(04) 






0.014 


m(fu) 




6.0 


6.0 



TABLE II: Total energies and magnetic moments of the 
atomic sites calculated for the A- and G-type antiferromag- 
netic states of the orthorhombic phases of Ca3Mn207. The 
energies are given in units of eV per half the FU (i.e., per 
Mn) with respect to those of the FM state. The magnetic 
moments are given in units of Bohr magneton (/is) 



Parameters 



A-type AF G-type AF 



AEaf~fm (eV) 


- 0.0082 


-0.0378 


fi{Mn) 


±2.613 


±2.678 


m(oi) 


±0.058 


~ 0.0 


m(02) 


0.0 


±0.03 


m(03) 


0.0 


±0.05 


m(04) 


0.0 


0.02 


m(fu) 


±5.168 


±5.380 




FIG. 1: (color online) Structures of (a) the minority tetrago- 
nal and (b) the majority orthorhombic phases of Ca3Mn207. 
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FIG. 2: (color online) Projected DOS plots of the atomic sites 
calculated for the NM states of (a) the minority tetragonal 
and (b) the majority orthorhombic phases of Ca3Mn207. 
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FIG. 3: (color online) Projected DOS plots of the Mn-3d or- 
bitals calculated for the NM states of (a) the minority tetrag- 
onal and (b) the majority orthorhombic phases of Ca3Mn207. 
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FIG. 4: (color online) ECOV plots of the Mn-0 bonds cal- 
culated for the NM states of (a) the minority tetragonal and 
(b) the majority orthorhombic phases of Ca3Mn207. 



(b) 



FIG. 5: (color online) Schematic representations of the Mn 
spin arrangements in the (a) FM, (b) the G-type AF, and (c) 
A- type AF states of Ca3Mn2 07. Blue and white spheres refer 
to up- and down- spin aUgnments respectively. 
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FIG. 6: (color online) FM Ca3Mn207. Site and spin pro- 
jected DOS (a) and chemical bonding (ECOV) for Mn-O in- 
teraction regrouping all four oxygen sublattices (b). Energy 
reference is with respect to the VB top Ev due to the closely 
insulating character. 
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FIG. 7: (color online) AF Ca3Mn207 : Site and spin pro- 
jected DOS (a) and band structure (b). Plots are for one 
magnetic sublattice. O-PDOS regroup all oxygen sites. 



